Ground state and constrained domain walls in Gd/Fe multilayers 
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The magnetic ground state of antiferromagnetically coupled Gd/Fe multilayers and the evolution 
of in-plane domain walls is modelled with micromagnetics. The twisted state is characterised by a 
rapid decrease of the interface angle with increasing magnetic field. We found that for certain ratios 
M Fe : M Gd , the twisted state is already present at low fields. However, the magnetic ground state 
is not only determined by the ratio M Fe : M Gd but also by the thicknesses of the layers, that is 
the total moments of the layer. The dependence of the magnetic ground state is explained by the 
amount of overlap of the domain walls at the interface. Thicker layers suppress the Fe aligned and 
the Gd aligned state in favour of the twisted state. Whereas ultrathin layers exclude the twisted 
state, since wider domain walls can not form in these ultrathin layers. 

PACS numbers: 



The study of spin-polarised transport is key for the un- 
derstanding of processes involved in spin-electronics. A 
particular interesting point is the effect of domain walls 
(DWs) on the transport properties of ferromagnetic ma- 
terials, as has been shown recently. 0, 0] This interest 
is based on the possibility of creating sharp boundaries 
(high angular rotation of the magnetisation in only few 
atomic layers) without any non-magnetic spacer. 

DWs have been studied in several systems. Domain 
walls can be measured in single crystal highly anisotropic 
thin films of SrRu03.[j| Or between features in a pat- 
terned Lao.7Cao.3Mn03 thin fihn.0 Another form of 
magnetic boundary has been observed in single crys- 
talline thin films of Fe3C>4, where at stacking faults the 
coupling is antiferromagnetic and 180° magnetic bound- 
aries are formed. Q It has been shown that domain walls 
can be formed in Gd/FM multilayer devices, where FM 
is a standard ferromagnetic metal like iron, cobalt or 
permalloy. 0, Q Recently, magnetoresistance measure- 
ments have been reported with current perpendicular to 
the in- plane DWs. 0,0 

From the point of view of transport properties of FM 
metals, Gd/FM multilayers are of special interest for 
three main reasons: firstly the Gd/FM interface is an 
antiferromagnetic (AFM) sharp boundary, so it is ideal 
for the study of spin diffusion and spin-flip. Secondly, 
there is a large resistance mismatch across the interface, 
making the system similar to metal-semiconductor inter- 
faces. Finally, in certain conditions as explained below, 
Gd/FM multilayers can support in-plane domain walls 
(DWs) close to the interface, so their transport proper- 
ties can be studied in detail. 

Current perpendicular-to-plane magnetoresistance ex- 
periments in Gd/Fe multilayers have shown that the re- 
sistance decreases with increased applied field due to a 
decrease in the interface angle c^.Q This is in qualitative 
agreement with magnetoresistance experiments in Fe304, 
where the resistance is a function of on controlled by the 
applied external field 5]. However, one of the samples in 



Ref. [21 showed an upturn in the resistance at large fields 
and it was speculated that this is related to the formation 
of DWs. 

The ground state domain walls at the Gd/Fe bound- 
ary have been magnetically characterised in bilayersj^ 
and multilayers too, but normally in quite thin layers. 
In these situations the interface and the AFM coupling 
play the most important roles in determining the mag- 
netic structure. Traditionally, the Fe aligned state is 
found for M Fe > M Gd and the Gd aligned state when 
M Fe < M Gd . Upon application of an external field H 
the twisted state with DWs at the interfaces is observed. 
Note the distinction between M Gd and M Gd . M Gd is 
the saturation magnetisation per volume and M Gd is the 
total moment of all Gd layers together which is thickness 
dependent. Most articles focus on the effect of changing 
ratio M Gd : M Fe with temperature, and therefore chang- 
ing ratio M Fe : M Gd . Since the Curie temperature of Gd 
is much smaller than that of Fe, M Gd decreases rapidly 
with respect to M Fe with increasing temperature. 

Micromagnetic simulation has been done to gain in- 
sight in the changes in the magnetic state with applied 
field. We study the effect of the thicknesses tFe and tad 
on the magnetic behaviour of Gd/Fe multilayers. We will 
show that there are distinct differences between samples 
with M Fe > M Gd and samples with M Fe « M Gd . We 
will also show that at constant ratio M Fe : M Gd , and 
therefore the magnetic ground state is significantly in- 
fluenced by the actual value of tad and tp e . We will 
present the phase diagram of the magnetic ground state 
as a function of ta and tp e . In ultra thin films the DWs 
on both side of the layer are not independent and interact 
via their surface tension. The Fe and Gd aligned state 
are completely suppressed for multilayers with tad > 16 
nm and tFe > 9 nm. 

In this paper we model Gd/Fe multilayers using the 
LLG micromagnetics software. The modelled multi- 
layers consist of three Fe layers, thickness tFe, and two 
Gd layers with thickness tad- The simulation volume is 
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10 x 10 x t t ot nm 3 , where ttot is the total thickness of 
the multilayer given by t to t = 3 tp e + 2 tad- The cell 
size in the simulation is 1 x 1 x 0.25 nm 3 , ensuring that 
the typical DWs found in this work are at least five to 
ten cells thick. In-plane magnetisation is obtained by 
periodic boundaries on the xz and yz planes. The mate- 
rials parameters for Fe and Gd that have been used are 
listed in Table HI effectively the simulations were done 
at T = K. The value used for M Gd is 25% smaller 
than the bulk value as reported in many experimental 
works. 0, 0, We will discuss the effect of reduced 
M Gd later on. The AFM Fe-Gd interface interaction is 
chosen to be Ai nt = —1.0 //erg cm -1 , between the ex- 
change coupling strengths of Fe and Gd. Furthermore, 
comparing work on the temperature dependence of the 
Co-Gd coupling and the Fe-Gd coupling indicates that 
Aint = —1-0 //erg cm -1 is a good estimate. 
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FIG. 1: Calculation of the DW width from the calculated 
magnetisation angle <j) versus layer depth t data. The drawn 
line is a fit to the data (o). The dashed line is the tangent 
at t' — 0, see text for explanation. The horizontal arrow 
indicates the DW width. 



TABLE I: Material parameters for Fe and Gd. 

Fe Gd 
M s (emu cm" 3 ) 1700 1508 
A (fierg cm" 1 ) 2.05 0.75 
K (//erg cm" 1 ) 1.4xl0 5 10 2 

(cubic) (uniaxial) 

For each simulation, a part of the hysteresis loop has 
been modelled from /j, H = 7.5 T to fi H = T. This 
prevented the problem of having to choose the correct 
initial state at /iqH = T. The initial magnetisation 
was aligned with the x-axis and the Fe (Gd) magneti- 
sation was parallel (antiparallel) with the positive di- 
rection. The field was applied in the positive y direc- 
tion. The initial field ^lqH — 7.5 T rotates the bulk of 
the layers parallel with the field; the Fe (Gd) boundary 
magnetisation however remains a component parallel to 
the positive (negative) x axis. All magnetisation vectors 
within a single xy simulation layer are identical and the 
z component is always negligible (< 10~ 4 Af s ). The angle 
between the magnetisation at layer depth t and the ap- 
plied field direction is <fr(t). Thus <j> = 0° corresponds to 
the magnetisation aligned with the y-axis and <f> = ±90° 
is aligned with the ir-axis. 

At low temperatures and zero field, each layer is ex- 
pected to be fully FM aligned, and antiparallel with re- 
spect to the adjacent layers, creating an interface angle 
a>i = 180°, with on = \4* Gd \ + \<j> F e |, where 4>b is the angle 
of the moment at the boundary layer. 

The in-plane DW widths at the Gd/Fe interfaces have 
been calculated as shown in Fig. ^ The orientation of 
the moments as a function of the distance to the interface 
is fitted to an exponential function of the form 4>{t') = 
aexp(6t'), where t' is the distance to the interface, a = 
4>{t' = 0) and is a fit parameter. The DW width 5 is 
the intersection of the tangent of this function at t' = 
with the = 0° axis. 

Three magnetic states have been identified in Gd/Fe 



multilayers. [tI lil HH IT3| Upon application of a small ex- 
ternal field H either the Fe layer or the Gd layer aligns 
with H . The ratio between the total magnetic moment 
per element M Gd : M Fe determines whether the Fe 
aligned (F) or Gd aligned (G) states prevails. For in- 
stance, Hosoito et al. report for a [Fe/Gd]i5 multilayer 
film that the G state is present for low temperatures 
T < 120 K where M Gd > M Fe , and the Fe aligned state 
prevails for high temperatures T >= 140 K, where M Gd 
has become smaller than M Fe since M Gd is significantly 
reduced. Both states are characterised by a* ~ 180° 
in small fields. A typical example of a micromagnetics 
simulation of the Fe aligned state is given in Fig. [3 show- 
ing a sample with tp e =10 nm and ta = 5 nm. This 
simulation will be referred to as the F state simulation 
in the remainder of the article. 
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FIG. 2: Angle <f> against the layer depth t for the simulation 
with tFe = 10 nm and tad = 5 nm. This is an example of 
the Fe aligned state. For small H the magnetisation in the Fe 
layer is fully aligned with the applied external field. At large 
H twisting occurs in the layers and a, is reduced. 
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Upon application of a large external field, most of each 
layer will align parallel with the field. Nevertheless, cj> will 
be nonzero near the interfaces. This is called the twisted 
(T) state. Fig. [21 sketches the alignment of the layers 
with the applied field (parallel to the y-axis) and how 
the magnetisation rotates away in the xy-plane near the 
interfaces to accommodate the AFM coupling, thereby 
creating in-plane DWs. 





H//y 



FIG. 3: Three dimensional view of the fioH = 7.5 T data of 
Fig-Hi F° r clarity only the data are plotted for 12.5 < t < 20 
nm, showing half of a Gd layer (thick arrows) and half of an 
Fe layer (thin). The magnetic field is parallel to the y-axis. 
At the interface, the Fe layer is -45° off the {/-axis and the Gd 
layer is +100° off the j/-axis. 

The twisted state originates from the balance between 
the AFM coupling at the interface (that tries to align one 
of the layers antiparallel to H) and the Zeeman energy 
(which tries to align all layers parallel to H). Fig.^Jshows 
the twisted state simulation in a sample with tp e =15 
nm and ted — 20 nm. The twisting of the FM layers 
is clearly present at all H > 0. Note that the layers at 
the edge of the samples rotate faster, which is due to 
the free surface. [l5j Also the high field data in Fig. 
shows the twisted phase. The lower exchange coupling of 
Gd allows a larger angle between neighbouring spins, i.e. 
favouring narrower DWs and DWs incorporating a larger 
twist, compared with the Fe layer. 

In experimental reports on Gd/Fe multilayers, the sat- 
uration magnetisation of the thin Gd layers was found to 
be reduced by 20% to 40%. [H GH E3 

We have studied 

the effect of reducing Mf d by comparing simulations of 
a tp e = tad = 30 nm multilayer. Two simulations have 
been performed, one with Mf d = 2010 emu cm -3 cor- 
responding to the bulk value and the other with Mf d 
= 1508 emu cm -3 in the range of the literature values 
found for thin Gd/Fe multilayers. The angles at the in- 
terface (fit, for both Gd and Fe layers, are extracted from 
the simulation and plotted in Fig. [SI 

Changing Mf d to 75% of the bulk value has only minor 
effects on the simulation. Experimental results indicate 
that the saturation magnetisation in thin Gd films might 
be reduced by 20% to 40%. Therefore M s Gd was reduced 
to 75% of the bulk value in all other simulations. 
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FIG. 4: 4> against t for the simulation with tFe =15 nm and 
tad = 20 nm. This simulation illustrates the twisted ground 
state. Even at small H on decreases rapidly with H and DWs 
are formed in both Fe and Gd layers. <j> c = 0° at relatively 
small H but 4>b 2> 0° for all H. Note that the twisting is more 
pronounced at low H in the top and bottom layer as recently 
reported by Haskel et al. 0] 
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FIG. 5: Angle at the boundary <j>b as a function of external 
field. Open symbols correspond to <j>^ d and solid symbols to 
Circles and triangles correspond to simulations using 
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1508 emu cm" 5 and Af s ° ' = 2010 emu cm , respec- 



We have also simulated a multilayer with a reduced in- 
terface coupling Ai n t — —Acd = —0.75 /xerg cm . The 
main features, relative to the "standard" simulation with 
Aint = —1-0 /icrg cm -1 , were strongly reduced ai and 
boundary angles but slightly higher DW widths for all 
fields. However the DW width actually increased slightly 
with increasing field at high fields. This is related to the 
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drastic decrease of a, and the very small twists in the 
DWs at those fields. 

Generally speaking the stronger the AFM coupling, the 
larger the interface angle and boundary angles. There- 
fore, the DWs will become narrower and have a larger 
twist. The narrower DWs will have less overlap and the 
twisted phase will still be present at smaller layer thick- 
nesses. 

We now look at the effect of sweeping the magnetic 
field on the magnetic state of the Gd/Fe multilayers. We 
have plotted the angles <p c , on, 4>b — <j>c and DW width 5 
as a function of H in Fig. El data is taken from the Fe 
aligned state simulation (Fig. and twisted state simu- 
lation (Fig. 2}. The main difference between the aligned 
states and the twisted state can be seen in the top-right 
panel of Fig. HJ In the twisted state simulation a.i de- 
creases rapidly with H , whereas in the Fe aligned state 
simulation on w 180° for small fields. The reduction of oti 
with increasing field is in agreement with the experimen- 
tal results reported on anti phase boundaries in F23O4 
and AFM coupled interfaces in Gd/Fe multilayers .Q 
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FIG. 6: Data extracted of the twisted state simulation with 
tFe = 15 nm and ted = 20 nm and the Fe aligned state sim- 
ulation with tpe = 10 nm and tad = 5 nm. The angle of 
the magnetic moment at the centre <f> c , the interface angle oti, 
the twist in the DW 4>b — 4>c and the width of the DW 8 are 
plotted vs the applied field H. The figure highlights the dif- 
ferent behaviour of the Fe and the Gd layer. Open circles (T 
- Gd) correspond to the Gd layer of the twisted state simula- 
tion; Closed circles (T - Fe) correspond to the Fe layer of the 
twisted state simulation. Open diamonds (F - Gd) represent 
the Gd layer of the Fe aligned state simulation and closed di- 
amonds (F - Fe) represent the Fe layer of the Fe aligned state 
simulation. 

The decrease of on in the twisted state simulation is 



accompanied by the formation of DWs at the interfaces 
of the thick layers. In the Fe aligned state simulation a*, 
4>b and 4> c are only weakly dependent on the applied mag- 
netic field. Each Gd or Fe layer remains almost perfectly 
aligned, i.e. no twisting in the layer. At about /-J.qH = 0.7 
T, there is a change in the gradient of on, 4>b and 4> c with 
H . The layers also start to twist and form wide DWs. 

In both simulations we find that with increasing field 
the DWs become narrower and on decreases. The curves 
for 5ad converge at fi H ~ 5 T. Apparently, the Gd DWs 
in both simulations become independent of ted at these 
high fields. The data for 8p e suggest that the same will 
happen in the Fe layers at even higher fields (noH > 8 
T). 

The samples in Fig. [f)| have been simulated whilst 
sweeping the field from /j,qH = 7.5 T to hqH = T. 
In Fig. [7| also the increasing field sweep is shown for the 
Fe aligned state simulation with tp e — 10 nm and tad 
= 5 nm, taking the magnetic configuration at H = as 
starting configuration. 
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FIG. 7: Data extracted of the F aligned simulation with tFe 
— 10 nm and ted = 5 nm. Data has been taken in decreas- 
ing field, followed by increasing field. Note the hysteretic 
behaviour at small fields. Y-axis symbols are the same as in 
Fig.© 

Looking at the behaviour at small fields, during the 
decreasing field sweep, the Fe (Gd) layers rotate slowly 
but steadily to be parallel (antiparallel) with the field di- 
rection at fi n H = T. In the increasing field sweep the 
layers initially stay parallel just as they are at /iqH = 
T. Between fi H = 2 T and fi Q H = 2.5 T, the mag- 
netic configuration abruptly changed to coincide with the 
decreasing field curve. Clearly the increasing field sim- 
ulation got trapped in a metastable state, whereas the 
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decreasing field sweep changed smoothly with H. 

The simulation for the twisted state sample (with 
tp e =15 nm and ted = 20 nm) shows no signif- 
icant difference between the decreasing and increasing 
field sweeps. In a second twisted state simulation, after 
the H || y decreasing field sweep, the increasing field is 
parallel to x. H is then parallel with either the Fe or the 
Gd magnetisation. During the increasing field sweep a 
similar hysteresis like transition is observed as in the Fe 
aligned state simulation. This demonstrates the impor- 
tance of the initial magnetic state and the direction of 
the field sweep. 

We have shown the typical simulations for the twisted 
state (tp e = 15 nm and tad = 20 nm) and the Fe aligned 
state {tp e = 10 nm and tod — 5 nm). We now investigate 
the effect of changing the layer thicknesses whilst keeping 
the ratio M Fe : M Gd constant. If the magnetic config- 
uration is only controlled by this ratio, then these simu- 
lations should exhibit qualitatively similar behaviour. In 
Fig. |S1 we plotted as a function of fioH the data extracted 
for different values of tp e = ted- 
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FIG. 8: The simulation of Gd/Fe multilayers for various ip e 
= tod- The simulations for tp e — tea, = 20 nm and tFe = ted 
= 30 nm are identical above fioH = 2 T. For clarity only the 
the data for the Gd layers are shown. Y-axis symbols are the 
same as in Fig. [S] 

At low fields and tp e — tod < 10 nm the F phase 
prevails, as expected since 

M Fe > M Gd j n these 

— tod multilayers. The prevalence of the F phase can 
be deduced from the small value of (f>b — 4> c , the slow 
decrease of «i, and the large value of (j) Gd ~^> 90°. Only 
for higher fields a transition to the twisted phase takes 
place and 014 starts to decrease fast. These simulations 
with tp e — tod < 10 nm are similar to the Fe aligned 



state simulation shown in Fig. [5] However, for larger 
tFe — tod the multilayer adopts the twisted phase at all 
fields even though M Fe > M Gd . As in the twisted state 
simulation (tp e = 15 nm and tod = 20 nm) on is much 
smaller than 180° at any fi^H > T and 3> <fi c for 
small fields. 

The different behaviour for small and large tp e = tod 
can be explained by the overlap of DWs in the centre 
of each layer. For high fields 25od « tod and the two 
DWs do not meet at the centre of the layer; they act 
independently. However at smaller fields (larger 8cd) and 
smaller ted the DWs do overlap. For instance, Fig. [S] 
shows that at [i H = 2.5 T the Gd DW is 4.8 nm wide for 
the multilayer with tp e — ted = 5 nm. The overlap ping 
DWs reduce the twist by means of the surface tension [l6| 
in the DW and the thin layers each become FM aligned. 
In Fig. El this scenario has been sketched. Since M Fe > 
M Gd the Fe layer will align parallel to the applied field, 
and (f> Gd will approach 180° at H = 0, resulting in the F 
phase. For thicker layers, the twist is removed, i.e. the 
layers become totally FM aligned, only when the field 
is very small during a decreasing magnetic field sweep. 
These small fields are insufficient to rotate the Fe layers 
parallel with the applied field. Consequently, both layers 
are perpendicular to the applied field direction, cj> c ~ 
±90° at H = and the twisted phase is the ground state 
for thicker layers. 
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FIG. 9: Sketch of the DW overlap. The DW width is reduced 
from the top to the bottom panel. The surface tension of the 
overlapping DWs enhances the reduction of the twist. 

The DW thickness depends on the layer thickness as 
can be observed in the bottom right corner of Fig. |SJ Sod 
increases with decreasing layer thickness at constant field. 
We also note that with decreasing layer thicknesses, the 
H—Scd curve converges at larger fields. This confirms the 
independent behaviour of narrow DWs in thicker layers 
as discussed above. 
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We have shown that there is a transition from the F 
phase at tp e = tod — 5 nm to the T phase at tp e = tad 
= 30 nm, even though the ratio M Fe : M Gd remained 
constant. Therefore, the magnetic ground state does not 
only depend on the ratio between the total magnetisa- 
tion but also on the individual layer thicknesses. We 
have investigated this dependency for a range of ^Fe's 
and tGd&- We have chosen to use a single variable as 
indicator for the magnetic ground state: 4> Gd at H = 0. 
If <j) Gd is between 45° and 135° the ground state of the 
simulation is designated as the T phase. Outside these 
boundaries either the F phase (<f) Gd > 135°) or the G 
phase {(f> Gd < 45°) is assigned. The phase diagram of 
the magnetic ground state as a function of tp e and tad 
is plotted in Fig. [TDJ 
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FIG. 10: Phase diagram as a function of the total Fe mag- 
netisation (bottom axis) and the total Gd magnetisation (left 
axis). The drawn lines indicate the boundaries between the 
Gd aligned (open squares), the twisted (grey triangles) and 
the Fe aligned (grey circles) states. The corresponding layer 
thicknesses are given on the top (tFe) and right (tad) axes. 

The phase diagram confirms that the boundaries be- 
tween the three magnetic ground states are not at fixed 
ratios of M Fe : M Gd . The data strongly suggests that 
the G phase is absent if tp e > 9 nm. One expects the G 
phase for multilayers with M Gd 3> M Fe but the phase 
diagram indicates that this is only the case if tFe < 9 nm. 
Above this value, no matter how thick the Gd layer is, 
the simulation will always yield the T phase. Similarly, 
the F phase appears to be absent if tad exceeds w 16 nm. 
We can conclude from these observations that when both 
layer thicknesses exceed these critical values only the T 
phase is present, independent on the ratio M^ e : M Gd . 

Below t to t = 20 nm, the fraction in the T phase drops 
rapidly with decreasing t to t- Already at high fields, the 
DWs overlap and reduce their twist. These thin layers 
do not allow the existence of DWs at small fields. This 
suppresses the T phase for ultra thin multilayers. 

In conclusion, using micromagnetics Gd/Fe multilay- 



ers consisting of three Fe layers and two Gd layers have 
been modelled to study the effect of tFe and tad on the 
magnetic ground state. The twisted phase can be dis- 
tinguished from the Fe and Gd aligned phases by the 
rapid decrease of on at small fields. At the same time 
wide DWs are formed at the edges of the thick layers. In 
both aligned states oti remains almost 180° in small fields. 
Each Gd or Fe layer remains almost perfectly aligned, 
i.e. no twisting in the layer. Above a threshold field, the 
aligned phases are transformed into the twisted phase. 
In all cases, at large fields the exchange constants keep 
Sfc > $Gd, with the twist in the Gd DW larger than in 
the Fe DW. 

When the magnetic field is applied parallel with the 
magnetisation of one of the layers in an increasing field 
sweep hysteretic behaviour can be observed. The sample 
will show a near perfect alignment of that layer with the 
field and a sharp transition to the twisted state is seen. 

We have found that the magnetic ground state depends 
both on the ratio between the total magnetisation for Fe 
and Gd and on the individual layer thicknesses. The 
different behaviour for small and large layer thicknesses 
can be explained by the amount of overlap of the DWs. 

The phase diagram shows that independent on the ra- 
tio M Fe : M Gd only the twisted phase is allowed when 
both layer thicknesses exceed their respective critical val- 
ues. For thin layers all the magnetic ground states are 
allowed and the ground state found depends on the ratio 
M Fe : M Gd . Ultra thin layers do not allow the existences 
of DWs at small fields. This suppresses the twisted phase. 
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